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Abstract. The existence of the limiting pair correlation for angles between reciprocal geodesies 
on the modular surface is established. An explicit formula is provided, which captures geometric 
information about the length of reciprocal geodesies, as well as arithmetic information about the 
associated reciprocal classes of binary quadratic forms. One striking feature is the absence of a gap 
beyond zero in the limiting distribution, contrasting with the analog Euclidean situation. 



1. Introduction 

Let EI denote the upper half-plane and T = PSL2(Z) the modular group. Consider the modular 
surface X = r\H, and let II : H — > X be the natural projection. The angles on the upper half 
plane EI considered in this paper are the same as the angles on X between the closed geodesies 
passing through LI(i) , and the image of the imaginary axis. These geodesies were first introduced 
in connection with the associated "self-inverse classes" of binary quadratic forms in the classical 
work of Fricke and Klein [8j p. 164], and the primitive geodesies among them were studied recently 
and called reciprocal geodesies by Sarnak [22]. The aim of this paper is to establish the existence 
of the pair correlation measure of their angles and to explicitly express it. 

For g E T, denote by 9 g G [— ir, it] the angle between the vertical geodesic [i, 0] and the geodesic 
ray [i,gi]. For zi,Z2 £ H, let d(z\,Z2) denote the hyperbolic distance, and set 

II^H 2 = 2coshd(i,gi) = a 2 + b 2 + c 2 + d 2 , for g= ( a c b d ) € SL 2 (M). 

It was proved by Nicholls [16] (see also [T71 Theorem 10.7.6]) that for any discrete subgroup T of 
finite covolume in PSL2(M), the angles # 7 are uniformly distributed, in the sense that for any fixed 
interval / C [— ir, n], 

#{ 7 er:e 7 €J,d(i,7i)<12} |/| 

Inn — p ; = — . 

r^oo #{ 7 e r : d(i, t») s: R} 2vr 

Effective estimates for the rate of convergence that allow one to take |/| X e~ cR as R — > oo for 
some constant c = cp > were proved for T = T(N) by one of us [lj, and in general situations 
by Risager and Truelsen [20J and by Gorodnik and Nevo [S]. Other related results concerning the 
uniform distribution of real parts of orbits in hyperbolic spaces were proved by Good [10], and 
more recently by Risager and Rudnick [19] . 

The statistics of spacings, such as the pair correlation or the nearest neighbor distribution (also 
known as the gap distribution) measure the fine structure of sequences of real numbers in a more 
subtle way than the classical Weyl uniform distribution. Very little is known about the spacing 
statistics of closed geodesies. In fact, the only result that we are aware of, due to Pollicott and 
Sharp |18j , concerns the correlation of differences of lengths of pairs of closed geodesies on a compact 
surface of negative curvature, ordered with respect to the word length on the fundamental group. 
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This paper investigates the pair correlation of angles # 7 with d(i, yi) ^ R, or equivalently with 



< Q 



2 cosh R as Q — > oo. As explained in Section 2l these are exactly the angles 



between reciprocal geodesies on the modular surface. 

The Euclidean analog of this problem considers the angles between the line segments connecting 
the origin (0,0) with the integral lattice points (m,n) with m 2 + n 2 ^ Q 2 as Q — > oo. The 
pair correlation of all lattice point angles (rays are counted with multiplicity) is plotted on the 
right in Figure [T] When only primitive lattice points are being considered (rays are counted with 
multiplicity one) , the problem reduces to the study of the sequence of Farey fractions with the L 2 
norm ||m/rt||| = m 2 + n 2 . Its pair correlation function is plotted on the left in Figure [l) When 
Farrey fractions are ordered by their denominator, the pair correlation is shown to exist and it is 
explicitly computed in [5]. A common important feature is the existence of a gap beyond zero for the 
pair correlation function. This is an ultimate reflection of the fact that the area of a nondegenerate 
triangle with integer vertices is at least 1/2. 

For the hyperbolic lattice centered at i, it is convenient to start with the (non-uniformly dis- 
tributed) numbers tan(# 7 /2) with multiplicities, rather than the angles # 7 themselves. Employing 
obvious symmetries explained in Section [3j it is further convenient to restrict to a set of repre- 
sentatives Ti consisting of matrices 7 such that the point 72 is in the first quadrant in Figure [2} 
The pair correlation measures of the finite set Tq of elements tan(# 7 /2) with 7 G and ||7|| Q 
(counted with multiplicities), respectively of the set 21q of elements 7 with 7 G Ti and ||7|| ^ Q 
(counted with multiplicities) are defined as 



#{(7,7 / )er i 2 :|| 7 ||,||y|Kg, 7 '/ 7 , ^ tan(0y/ 2 ) -tan(0 7 /2) < j^} 

Rq(0 = — 



4(0 



Br 



#{(7,Y) G T 2 : || 7 ||, IIVII ^,7^7,0^ £(0y - 7 ) ^ ^} 



Br 



where Bq ~ |Q 2 denotes the number of elements 7 G Tj with ||7|| Q. 




Figure 1. The pair correlation functions (left) and (right), plotted in grey, 
compared with the pair correlation function of Farey fractions with L 2 norm, (left) 
and of the angles (with multiplicities) of lattice points in Euclidean balls (right). 
The graphs are obtained by counting the pairs in their definition, using Q = 4000 
for which Bq = 6000203. We used Magma [H] for the numerical computations, and 
SAGE [21] for plotting the graphs. 
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One striking feature, shown by numerical calculations in Figure [TJ points to the absence of a gap 
beyond zero in the limiting distribution, in contrast with the Euclidean situation described above. 

The main result of this paper is the proof of existence and explicit computation of the pair 
correlation measures -Rf and i?f , given by 



2 ' 

flf(O = ^((0^]):=Jim Rl(0, /#(0=i#((<U]):-,jini /?*,(<)■ M.l) 

thus answering a question raised in pQ . 

To give a precise statement consider 6, the free semigroup on two generators L = (} 5) anc ^ 
R = ( J i )• Repeated application of the Euclidean algorithm shows that 6 U {/} coincides with the 
set of matrices in SL2(Z) with positive entries. The explicit formula for R\ (£) is given as a series of 
volumes summed over 6, plus a finite sum of volumes, and it is stated in Theorem [2] of Section 7. 
The formula for i?f (£) leads to an explicit formula for R^ (£), which we state here, partly because 
the pair correlation function for the angles # 7 is more interesting, being equidistributed, and partly 
because the formula we obtain is simpler. 

Theorem 1. The pair correlation measure R% exists on [0, oo), and it is given by the C l function 

^(i^) = 3S2)( E B ««>+ E E A «M))- d.2) 

v 7 sw \Mee £e[o,£/2)ire[i,£/2) 

For M e <S, letting U M = ||M|| 2 /\/ll M ll 4 ~ 4 > °m as above, and f + = 

n - f (l/Vl|M|| 4 -4-sin(20-^ / )/e 



4 7 Um + cos (26-6 m ) 




For integers K G [l,£/2), ^ € [0,£/2) we ftcwe 



= / V4 ^ (200k'*) colt' 







wni/i ^4i^(£,t) t/ie area of the region defined by, setting e = (cos 9, sin 0). 



{l,Zf(^)+2^ +1 (^)} 



max 



where Fxi(0) := cot6> + V- , 7 — , ,- fl . r , jfl , H 7 — — - r and the linear forms LAx,y) 

K ' (V > ^-i( e ' 9 )^( e8S ) L,(e^)(L|( e ^)+L| +1 (e^)) ^ U ' ^ 

are defined in terms of the triangle transformation T : [0, l] 2 -4 [0, l] 2 , T(x,y) = (y, [^~^]y — %) 
by iteration: T l = (_Lj_i,Lj) /or 1 < i < I and L^ + i = KL^ — 

Rates of convergence in (1.1) are effectively described in the proof of Theorem [2] and in Propo- 
sition [T4l 



When £ ^ 2, the second sum in (1.2) disappears and the derivative B' M {£) is explicitly computed 
16, yielding an explicit formula for the pair correlation density function g^H) = ~rtr^ 



in Lemma 



which matches the graph in Figure [TJ 
Corollary 1. For < £ ^ 2 we have 



iw ) 3 ? £?e \\\Mf + ^\\M\\*-4-e 



A formula valid for < £ ^ 4 is given in (8.16) after computing A' K0 (£). 
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X 



_3_ 



The computation is performed in £8.2 and it identifies the first spike in the graph of gf{x) at 



y/E. An explicit formula for the pair correlation density g%(x) valid for all x, and also when 
the point i is replaced by the other elliptic point p = e 7 ™/ 3 , will be given in [3] . 

Since the series in Corollary [l] is dominated by the absolutely convergent £ 2 ||M|| -4 , we can 
take the limit as £ — > 0: 

g f(0) = \ y ( l|M| ' 2 - l) = 0.7015... 




Remarkably the previous two formulas, as well as (1.2) for £ ^ 2, can be written geometrically as 
a sum over the primitive closed geodesies C on X which pass through the point II (i), where the 
summand depends only on the length i(C): 

92 (°) = 3 e n£(C) _ l ■ 

C n^l 

This is proved in Section [2j where we also give an arithmetic version based on an explicit description 
of the reciprocal geodesies C due to Sarnak [22] . 

For the rest of the introduction we sketch the main ideas behind the proof, describing also the 
organization of the article. After reducing to angles in the first quadrant in Section [3| we show 
that the pair correlation of the quantities = tan(# 7 /2) is identical to that of $(7) = Re(7i). 
We are led to estimating the cardinality 7Zq(£,) of the set 

{(7,70 G T : 2 : || 7 ||, IIVH < Q, 7 V 7, < $( 7 ') - $(7) < J-}- 



For 7 = I ^, with nonnegative entries, ||7|| < Q, and g, g > 0, consider the associated Farey 

interval \p/q,p' /</], which contains $(7). In Section |4j we break the set of pairs (7,7') above in 
two parts, depending on whether one of the associated Farey intervals contains the other, or the 
two intervals intersect at most at one endpoint. In the first case we have 7 = j'M or 7' = 7M with 
M £ &, while in the second we have a similar relation depending on the number £ of consecutive 
Farey fractions there are between the two intervals. The first case contributes to the series over (3 in 
( |1.2[ ), while the second case contributes to the sum over K, L The triangle map T, first introduced 
in [3], makes its appearance in the second case, being related to the denominator of the successor 
function for Farey fractions. 

To estimate the number of pairs (7, 7M) in the first case, a key observation is that for each 



MET there exists a function Em(x,]j), given by (5.1), such that 



*(7)-$( 7 M) =S M (q',q) 

for 7 as above. Together with estimates for the number of points in two dimensional regions based 
on bounds on Kloosterman sums (Lemma|7|, this allows us to estimate the number of pairs (7, jM) 



with fixed M € S, in terms of the volume of a three dimensional body Sm,§ given in (7.14). The 
absence of a gap beyond zero in the pair correlation measure arises as a result of this estimate. The 
details of the calculation are given in Section [TJ giving an explicit formula for i?2 (Theorem [2]). 

Finally in Section [8] we pass to the pair correlation of the angles # 7 , obtaining the formulas of 
Theorem [T] and Corollary [TJ 

2. Reciprocal geodesics on the modular surface 



In this section we recall the definition of reciprocal geodesics and explain how the pair correlation 
of the angles they make with the imaginary axis is related to the pair correlation considered in the 
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introduction. We also show that the sums over the semigroup <3 appearing in the introduction can 
be expressed geometrically in terms of sums over primitive reciprocal geodesies. A description of 
the trajectory of reciprocal geodesies on the fundamental domain seems to have first appeared in 
the classical work of Fricke and Klein [HJ p. 164], where it is shown that they consist of two closed 
loops, one the reverse of the other. There the terminology "sich selbst inverse Classe" is used for the 
equivalence classes of quadratic forms corresponding to reciprocal conjugacy classes of hyperbolic 
matrices. 

Oriented closed geodesies on X are in one to one correspondence with conjugacy classes {7} of 
hyperbolic elements 7 G T. To a hyperbolic element 7 G F one attaches its axis a 7 on H, namely 
the semicircle whose endpoints are the fixed points of 7 on the real axis. The part of the semicircle 
between zq and 7^0 , for any zo G a 7 , projects to a closed geodesic on X, with multiplicity one if 
only if 7 is a primitive matrix (not a power of another hyperbolic element of T). The group that 
fixes the semicircle a 7 (or equivalently its endpoints on the real axis) is generated by one primitive 
element 70. 

We are concerned with (oriented) closed geodesies passing through H(i) on X. Since the axis 
of a hyperbolic element A passes through i if and only if A is symmetric, the closed geodesies 
passing through H(i) correspond to the set 1Z of hyperbolic conjugacy classes {7} which contain a 
symmetric matrix. The latter are exactly the reciprocal geodesies considered by P. Sarnak in |22j, 
where only primitive geodesies are considered. 

The reciprocal geodesies can be parameterized in a 2-1 manner by the set 6cT, defined in the 
introduction, which consists of matrices distinct from the identity with nonnegative entries. To 
describe this correspondence, let A C T be the set of symmetric hyperbolic matrices with positive 
entries. Then we have maps 

e^A^n (2.1) 

where the first map takes 7 G & to A = 77*, and the second takes the hyperbolic symmetric A to 
its conjugacy class {A}. The first map is bijective, while the second is two-to-one and onto, as it 
follows from [22j. More precisely, if A = 77* G A is a primitive matrix, then B = 7*7 ^ A is the 
only other matrix in A conjugate with A, and {A n } = {B n } for all n ^ 0. 

Note also that ||7|| 2 = 1^(77*), and if A is hyperbolic with Tr(A) = T, then the length of the 
geodesic associated to {,4} is 2\nN(A) with N(A) = T+v/ J 2 ~ 4 . 

We need the following: 

Lemma 2. Let A G T be a hyperbolic symmetric matrix and let 7 G T such that A = 77*. Then 
the point *yi is halfway (in hyperbolic distance) between i and Ai on the axis of A. 

Proof. We have d{i,"fi) = d(i,7*i) = d("fi,Ai) where the first equality follows from the hyperbolic 



distance formula and the second since T acts by isometries on EI. Using formula (3.2), one checks 



that the angles of i, 72 and i, Ai are equal, hence 72 is indeed on the axis of A. □ 

We can now explain the connection between the angles # 7 in the first and second quadrant in 
Figure [2J and the angles made by the reciprocal geodesies with the image TL(i — > ioo) = H(i — > 0). 
Namely, points in the first and second quadrant are parameterized by "fi with 7 G 6, and by the 
lemma the reciprocal geodesic corresponding to A = 77* G A consists of the loop LT(i — > ji), 
followed by TL(i — > 7*i) (which is the same as the reverse of the first loop). Therefore to each 
reciprocal geodesic corresponding to A = 77* G A correspond two angles, those attached to 7^ and 
7*z in Figure [2J measured in the first or second quadrant so that all angles are between and ir/2. 

In conclusion the angles made by the reciprocal geodesies on X with the fixed direction Tl(i —} 
ioo) consist of the angles in the first quadrant considered before, each appearing twice. Ordering 
the points "fi in the first quadrant by ||7|| corresponds to ordering the geodesies by their length. 
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Therefore the pair correlation measure of the angles of reciprocal geodesies is 2^(^/2), where K\ 
was defined in the introduction. 



The parametrization (2.1) of reciprocal geodesies allows one to rewrite the series appearing in 
the formula for <7<f(0) m the introduction, as a series over the primitive reciprocal classes lZ pnm : 

^ ( JtMW^Z ~ 1 ) = ^ N(A) 2 - 1 = 4 ^ ^ NH) 2n - 1 ' 
A/e6\vll iW ll * / AeA y ' { 7 }e7ep™n>i w/ 

where we have used the fact that for a hyperbolic matrix A of trace T we have 

VT 2 -4 = JV(A) - iV(A)- 1 and iV(A n ) = N{A) n . 

One can rewrite the sum further using the arithmetic description of primitive reciprocal geodesies 
given in |22j . Namely, let T>ti be the set of nonsquare positive discriminants 2 a D' with a G {0, 2, 3} 
and D' odd divisible only by primes p = 1 (mod 4). Then the set of primitive reciprocal classes 
7£ pnm decomposes as a disjoint union of finite sets: 

^prim = JJ ^prim 

with |7£^ nm | = z^(gQ, the number of genera of binary quadratic forms of discriminant d. For d G D^, 
z/(d) equals 2 A_1 , or respectively 2 A depending on whether 8 \ d, or respectively 8 | (i, and A is the 
number of distinct odd prime factors of d. Each class {7} G TZ^ nm has 

^(7) = "d = ^ 

with (iio,wo) the minimal positive solution to Pell's equation u 2 — dv 2 = 4. We then have 

2^ 2^ iV( 7 )2n _l 2^ a 2n _ y 



In the same way, by Lemma 12 the pair correlation measure -Rf (£) m Theorem 1 can be written 



for £ ^ 1 as a sum over classes {7} G 7£ pnm , where each summand depends only on £ and N(j). 

3. Reduction to the first quadrant 

In this section we establish notation in use throughout the paper, and we reduce the pair corre- 
lation problem to considering angles in the first quadrant in Figure [2} A similar reduction can be 
found in [6], in the context of visibility problems for the hyperbolic lattice centered at i. 

For each g = (" ^) G SL2 (M) define the quantities 

X g = a 2 + b 2 , Y g = c 2 + d 2 , Z g = ac + bd, T g = X g + Y g = \\g\\ 2 , 

<%) = Re( 9 i) = ^, e^^-^Tf^l). (3 ' 1} 

Note that X g Y g - Z 2 = 1. 

A direct computation yields that the angle 8 g G [— vr,7r] between the vertical geodesic [i,0] and 
the geodesic ray [i,gi] is given by 

* (9): . tm (|).i^_A_. ,3.2) 
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Figure 2. Two symmetric geodesies through i 



The upper half-plane HI is partitioned into the following four quadrants: 

I = {z G H : Rez > 0, \z\ < 1}, II = {z G M : Re z > 0, |z| > 1}, 
III = {z G H : Rez < 0, |z| > 1}, IV = {z G M : Rez < 0, |z| < 1}. 

We denote 7 = ("J) ,7 = (ba)' s= (So)- Let 7 € F, 7 ^ For 71 to be in the right 
half-plane we need Re(7«) > 0. This is equivalent with ac + bd > and implies ac ^ 0, bd ^ 
because abed = be + (6c) 2 ^ 0. Since ac ^ without loss of generality we will assume a ^ 
and c ^ (otherwise consider —7 instead). Without loss of generality assume b ^ 0, d ^ as 
well (otherwise can consider 7s = (Z^c) i ns tead since 'ji = •ysi), so we can assume 7 has only 
nonnegative entries. 

If a,b,c,d ^ and ad — be = 1, then - and I are both ^ 1 or both 1 (since open intervals 
between consecutive Farey fractions are either nonintersecting or one contains the other). Since 
yi G I <^=> a 2 + b 2 < c 2 + d 2 , it follows that both | and 3 are ^ 1 for 7? G I. We conclude that 
among the eight matrices ±7, ±75, ±7, ±7s, which have symmetric angles (see Figure [2]), the one 
for which 7* is in quadrant I can be chosen such that 



a, b, c, d ^ and 



b a 

< -= < - < 1 
a c 



The set of such matrices 7 is denoted Ti. 

Consider the subset £Rq of T\ consisting of matrices with entries at most Q, that is 



P P 



ET:0^p,p',q,q' ^ Q, 



p p 
< . 



K 1 



and its subset 9\q consisting of those 7 with ||7|| ^ Q. The cardinality Bq of %\q is estimated 
in Corollary 8] as Bq ~ in agreement with formula (58) in [22] for the number of reciprocal 
geodesies of length at most x = Q 2 . 

Let J-Q be the set of Farey fractions p/q with < p < q < Q and (p, q) = 1. The Farey tessellation 
consists of semicircles on the upper half plane connecting Farey fractions ^ p/q < p'/q 1 ^ 1 with 
p'q — pq' = 1 (see Figure [3]). We associate to matrices 7 G 9\q with entries as above the arc in the 
Farey tessellation connecting p/q and p'/q'. We clearly have 

#<R Q = 2#T Q - 3 = ^ + O(QlnQ). 
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4. The coincidence of the pair correlations of $ and ^ 

In this section we show that the limiting pair correlations of the sets {^(7)} and |<I>(7)} ordered 
by UtII ~~ * 00 coincide. The proof uses properties of the Farey tessellation, via the correspondence 
between elements of D\q and arcs in the Farey tessellation defined at the end of Section [3j 

For 7 = (J p ) G SH Q , set 7 _ = §, 7+ = J- From §U$, Q we have: 



Z n 1 



U) ™ YgtfYg-l) IMI 4 ' 1 ' 

7_ < $(7) < $(7) < 7+. (4.2) 
Denote by 7£q(£), respectively TZq(^), the number of pairs (7,7') G 9^ such that < #(7) - 
*(t0 ^ q2 5 respectively < $(7) - $(7') < gj. For fixed /3 G (§, l), consider also 

MQM ■= {(7,7) G : l*(7) " * (7)1 < ^, IMI < 
and the similarly defined Mq g p ■ The trivial inequality 

H%(t) < 2A/|,^ + #{(7, 7') G ^ : |*( 7 ) - *(Y)| < INI, IIVII > Q' 3 } 
and the estimate in (4.1) show that there exists a universal constant k > such that 

< + #{(7,7') G 9$ : 7 * Y, < *( 7 ) - *( 7 ') < ^ + J^}, 

showing that 

ttg(0 < 2N^ p + 1I%{2kQ^) + H*(Z + 2kQ 2 -^). (4.3) 
In a similar way we show that 

7^(0 < 2< fi/3 + 7e*(2 K Q 2 - 3 ' 3 ) +7 e*(e + 2 K Q 2 ^). (4.4) 

We first prove that Mq ^ p and Mq ^ p are much smaller than Q 2 . For this goal and for latter use, 
it is important to divide pairs (7, 7') G £Kq in three cases, depending on the relative position of their 
associated arcs in the Farey tessellation (it is well known that two arcs in the Farey tessellation are 
nonintersecting) : 

(i) The arcs corresponding to 7 and 7' are exterior, i.e. 7+ ^ 7L or 7^ ^ 7 _. 

(ii) i ~ 7, i-e. 7- < 7- < 7+ < 7+- 

(iii) 7 £ 7', i.e. 7L 5^ 7_ < 7+ ^ 7+ . 

Proposition 3. M^p < Q 1+/3 lnQ and M^p < Q 1+/3 lnQ. 

Proof. Mq £ p and A/q ^ ^ are increasing as an effect of enlarging !lKq to SHq, so for this proof we will 
replace 9\q by 9\q. We only consider Mq ^ p here. The proof for the bound on Mq t a is identical. 
Both rely on ( |4~lj ) and ( |42| ). 

Denote K =[£} + 1. Upon ( |42| ) and |r' - r| ^ Vr, r' G Jq, r / r', it follows that if 7+ < i_ 

and - <I>(-x) I < qz, then #(Jq D [7+, 7L]) < -fiT + 1. In particular 7^. = 7+ when < £ < 1. 

We now consider the three cases enumerated before the statement of the proposition. 

(i) The arcs corresponding to 7 and 7' are exterior. Without loss of generality assume 7+ ^ 7L • 

If i is such that 7+ = 7$, the i th element of J-q, then = ji +r = ^ 1±J - for some ^ r < K. 

p' p' 
The equality p + q_ — pLq+ = 1 shows that if j'_ = is fixed, then q' + (and therefore 7+ = ^) 
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is uniquely determined in intervals of length ^ q'_. Since q± ^ Q, it follows that the number of 
choices for q', is actually ^ S- + 1 = — h 1. 

When < ^ < 1 one must have ^y'_ = 7+. Knowing g_ and q+ would uniquely determine the 
matrix 7. Then there will be at most — + 1 choices for 7', so the total contribution of this case to 

< E E 

When £ ^ 1 denote by ft+i, . . . , the denominators of 7j, 7i+i 5 • • • , 7i+_ft"- As noticed in [XT] . 
= ~ 3r As in 13 consider nfay) := f 1 ^] and 7^ = {(x,y) G (0, l] 2 : x + y > 

= fc}. Let Q large enough so that So := Q 13 ^ 1 < 2 k+3 • Then % < So and it is plain 
(cf. also m that «Ji > 1 - S and «(§, ^i) = • • • = *(*^, = 2. Hence (|, Sgi) G 7i 

and (5i±i,^±2),..., (^t^, 211 ^) G 75, showing in particular that min{%+i, . . . , q i+K } > f- 
Therefore max ( -Q— , . . . , — ) < 3 and the contribution of this case to A/",? t « is 

< E E 4K«^. 



Q 



+ 1 ) <Q 1+/3 lnQ. 




Figure 3. The Farey tessellation. 

(ii) 7' ;$ 7. Let i be the unique index for which 7^ < $(7) < 7i+i with ji < successive 
elements in Tq. Since ($(7') — ^(7)! ^ qi, either < $(7) < 7+ or there exists ^ r ^ K with 
7^_ = 7j_ r or with 7^ = 7i+ r . In both situations the arc corresponding to the matrix 7' will cross at 
least one of the vertical lines above ji-K, ■ ■ ■ , liili+i-, ■ ■ ■ , li+K- A glance at the Farey tessellation 
provides an upper bound for this number N 7> k of arcs 7' G 9\q- Actually one sees that the set 
C-y l consisting of 2 + 2 2 + • • • + 2 L arcs obtained from 7 by iterating the mediant construction 
L=[ min{^,g + }] + * times (7 is not enclosed in C 7i l) contains the set {7' G 91q : 7' ^ 7,7' 7^ 7}. 
The former set contains at most L arcs that are intersected by each vertical direction, and so 
Nj,k ^ (2AT + 1)L. Therefore, the contribution of this case to A/q ^ p is (first choose 7, then 7') 

< E E ^^i-^]^)^^- 
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(iii) 7 5; 7'. We necessarily have 7 = 7'M, with M G G. In particular this yields 7j_ G -^q/s- 
Considering the sub-tessellation defined only by arcs connecting points from Tqp , one sees that the 

number of arcs intersected by a vertical line x = a with 7_ = |<a<7 + = | 7 ,7 = (t_,7 + ) G Tqh 
is equal to s(q, q'), the sum of digits in the continued fraction expansion of i < 1 when q < q' , and 
respectively to s(q',q) when q' < q. A result from [25j yields in particular that 

£ «Q 2/? In 2 Q, 

and therefore 

#{ (7, V) G : 7 £ YX 1 + 2 £ g') « g 2/3 In 2 Q. 

0<q<q%QP 

This completes the proof of the proposition. □ 



Proposition [3] and inequalities ( |4.3| ) and (4.4) provide 



Corollary 4. For each /3 G (1,1 



3> -v> 



^q(0 = ^q(^ + 0(Q 2 - 3/3 )) +^S(0(Q 2 ^)) +0(Q 1+/3 lnQ). 

5. A DECOMPOSITION OF THE PAIR CORRELATION OF {$(7)} 

To estimate IZq(^), recall the correspondence between elements of 5Hq and arcs in the Farey 
tessellation from the end of Section |3| We consider the following two possibilities for the arcs 
associated to a pair (7,7') G 9\q: 

(i) One of the arcs corresponding to 7 and 7' contains the other. 

(ii) The arcs corresponding to 7 and 7' are exterior (possibly tangent). 

Denoting by Rq(Q, -Rg n (£) the number of pairs in each case we have 

^5(0 = *5(e)+*8 n (fl- 

5.1. One of the arcs contains the other. In this case we have either 7 = 7'M or 7' = 7M with 
M G 6 (see also Figure [I]). For each M G T define 

„ / v ~ -X"m) + (^ 2 ~ ^ 2 )^M (t .-is 

*M[x,y) ix 2 + y 2 )ix 2 XM+ y2 YM + 2xyZ M y ^ 

where Xm,Ym, Zm are defined in ( |3.1[ ). A direct calculation leads for 7 = (^ P q > ^ to 

$( 7 )-$(7M) = SM(g',g) (5.2) 

proving the next statement. 

Lemma 5. The number of pairs (7,7') G 9tg, 7 / 7' , ^ $(7) — $(7') ^ qs an d 7 « 7' or 
7' k, 7 ^ e notation introduced before Prop. [5]) is given by 

RUO = #{(7,7M) G 5% : 7 = (j! , M G 6, \E M (q',q)\ < 
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Figure 4. The case 7' 



< 



5.2. Exterior arcs. In this case we have 7,7' S !lKq, 7L. #s 7+. Let £ > be the number of 
Farey arcs in Tq connecting the arcs corresponding to 7,7' (see Figure |5j). In other words, writing 

7 = (g'g)'7 /= f qe+i qe ) 1 we have that |a — 2_ ) Ei )---) £i. are consecutive elements in Tq . Setting 

also H=i := 2 it follows that q { = if^-i - %_ 2 , where #j e N, i = 1, . . . ,£, and *Q = \ Q t 9i ~ 2 ] for 
2 ^ » ^ 1 

The fractions ^i, are not necessarily consecutive in Tq, but we have qn + \ = — qe~i, 
K ^ K l+1 = [2±fc]. It follows that Y = 7M with M = ( K \ J ) • • • ( f< J ) ( * ^ 




FIGURE 5. The case when 7 and 7' are exterior 



We have t < £ because 



It is also plain to see that 



1 1 

SKY) - $(7) > £ - — ; > 



1=1 



<?i-i<7i ^ <2 2 ' 



<?' m sV + s' 2 )' m % %(«| + gl + i)' 

The last equality in (5.3) and g| + ^ Q 2 yield for f ^ 1 



X > <w 7 ') _ $ (7 ) > 1 + ««±i 

q 2 ^ 2 +<? 2 +i; 



while if ^ = we have 



> 1 ifrft - g-i = K | Q 2 - > ^ 

17 ' [V {q 2 + q' 2 ){q 12 + q 2 ) Q 2 ' 



(5.3) 
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showing that K < £. Notice also that (5.3) yields 

Let T = {(x, y) G (0, l] 2 : x + y > 1} and consider the map 
T:(0,1] 2 ^T, T(x,y)=(y, 



Qe+i 



2 + ll+iY 





"1 + x~ 








y-x\ 




y 





Its restriction to T is a area-preserving bijection of T |3j. Denote T' L = (Lj_i, Lj) and K{ = 2 
if £ = 0, 1, . . . , £, Kp + i = K, and L^ + i = KL^ — Lg_i, so one has 

< Li(x, y) ^ 1, i ^ 0, Li-t(x,y) + Li(x,y) > 1, i = l,...,£, 
L-i(x,y) = x, L (x,y) = y, 

Li(x,y) = K i (x,y)L i ^ 1 {x,y) - L^ 2 (x,y), i = 0, 1, ...,£+ 1, 



(q l -l,q i )=QT i 



Q_ Q_ 
Q'Q 



i = 0,l,...,i 



(5.4) 



q l+ i = Kqn - q^i = Q \^KL e (^, ^ ) - L^_i ^ 9 9 



Define also the function 

12 



T^:(0,l]^(0,oo), 



1 

Z72\ + X] 777 



+ 



(5.5) 



" L (L 2 _ 1 + Lq) ^ Li-iLi L e (Lj + L 2 +1 ) ' 

We proved the following statement. 

Lemma 6. The number i?Q n (£) of pairs (7, 7') of exterior (possibly tangent) arcs in IKq for which 
< $(t / ) - $(7) ^ ^2 «s given by 



O^p^q, O^p' ^q', p'q- pq' = 1 
p p\ P 2 + p' 2 + q 2 + q' 2 HQ 2 , 0<Kq t -q e -i^Q 
' q) : pj + q 2 e + (Kpe - Pe-1) 2 + (Kq e - %_i) 2 H Q 2 



(5.6) 



where the sums are over integers in the given intervals, and q-\ = q,qo = q' . 

6. A LATTICE POINT ESTIMATE 

Lemma 7. Suppose that is a region in M? of area A(£2) and rectifiable boundary of length £(d£l). 
For every integer r with (r, q) = 1 and 1 ^ L H q 



Mi 



(), /iP := #{(a,b) £ nnZ 2 : ab = r (mod q)} - 



A(n)+£ n ,L 



where, for each e > 0, 
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Proof. Replacing Z 2 by LZ 2 in the estimate (for a proof see [151 Thm. 5.9]) 

{(m,n) G Z 2 : (m,m + 1) x (n,n + l)fl9fl} <C l+*(9fi), 

we find that the number of squares S m<n = [Lm, L(m + 1)] x [Ln, L(n + 1)] with S m:n D dQ, ^ is 
<C l + ^£(dn). Therefore 

#{(m, n) G Z 2 : (Lm, L(m + 1)) X (Ln, L(n + 1)) C ft} = ^ + O (l + . 

Weil's estimates on Kloosterman sums |23| extended to composite moduli in [12] and [7] show that 
each such square contains ^r-L 2 + E (q 1 ^ 2+£ ) pairs of integers (a, 6) with ab = r (mod q) (see, e.g. 
[21 Lemma 1.7] for details). Combining these two estimates we find 

«w = (4*> + o(i + + otf/-,) = *>, (si) + ^, 

as desired. □ 
Corollary 8. = ^ + £ (Q 11/6+e ). 

Proof. Note first that one can substitute ~ for p' = 1+ J q in the definition of 91q, replacing the 
inequality ||7|| 2 ^ Q 2 by (q 2 + q' 2 )(q 2 + p 2 ) ^ Q 2 g 2 without altering the error term. Applying 
Lemmamto VL q = {(u,v) G [0, g] x [0, Q] : (g 2 + u 2 )(g 2 + v 2 ) ^ Q 2 q 2 } and L = g 5 / 6 we infer 

#^ = E^ M -^ 1 +o £ (Q 11/6+e ). 

^ g g y 

Standard Mobius summation (see, e.g., J2J Lemma 2.3]) applied to the function h(q) = ^A(Q q ) 
with ||/i||oo ^ Q and the change of variable (q',u,v) = (Qx,Qxy,Qz) further yield 

#^Q = ^yVol(5) + £ (g 11 / 6+e ), 

where 

S = {(x, y, z) G [0, l] 3 : (1 + y 2 )(x 2 + z 2 ) < 1}. 
The substitution y = tan yields 

[■it/A in 2 

Vol(S) = / -^^({(^)e[0,l] 2 :* 2 + ^ cos 2 #}) = -, 
Jq COS U id 

completing the proof of the corollary. □ 

The error bound in Corollary [8] can be improved using spectral methods (see Corollary 12.2 in 
Iwaniec's book [13]). We have given the proof since it is the prototype of applying Lemma [7] to the 
counting problems of the next section. 

7. Pair correlation of {$(7)} 

The main result of this section is Theorem [2j where we obtain an explicit formula for the pair 
correlation of the quantities {$(7)} in terms of volumes of three dimensional bodies. The discussion 
is divided in two cases, as in Section [5j 
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7.1. One of the arcs contains the other. The formula for Rq in Lemma [5] provides 

Mee 

where Mm,q(0 denotes the number of matrices 7 = ( ^, ^ for which 

O^p^g, 0<p'W, p'q-pq' = l, \E M (q',q)\^-^, WKQ. (7.2) 
The first goal is to replace in (7.2) the inequality ||7M|| ^ Q by a more tractable one. Taking 

/2 



7 = (yl and substituting p = p q , we write, using the notation (3.1 ): 



|| 7 M||' 



p_ 

n>1 



+ l) W 2 X M + q*Y M + 2qq'Z M ) + ^ - 2 ^+/^ . (7.3) 

Let A/j\,/(Q,0 denote the number of integer triples (q',q,p') such that 
sC p' sC g' sC Q, 0<q^Q, p'q=l (mod g'), 

QAf q)\ < ^2, /^M + ^Ym + 2gg'Z M «S ^l q ,2 - ^'^ 

If Ya/ ^ Q 2c ° for some constant Co G (|, l), then (7.3) yields 

■A/at,q(£) = A/m (Qa + OCQ 00-1 ))^) • (7-5) 
Next we show that Nm,q{£) = when max{JM, Ym} ^ Q 2c ° and Q is large enough. 

Lemma 9. Let cq G (|>l)- There exists Qo(Q such that N"m,q(0 = for every M G 6 
max{X M ,y M } ^ Q 2c « and Q ^ Q (0- 

Proof. We show there are no coprime positive integer lattice points (q',q) for which 

|Sm(</, </)l ^ 7%, Y 7 m = 9 /2 X M + g 2 Y M + 2qq'Z M < Q 2 . (7.6) 



Suppose (g',g) is such a point, write in polar coordinates q'i + q = (q,q') = (r cos 9 , r sin 9) , 9 G 
(0,f), and consider (X,Y,Z) = (X M ,Y M ,Z M ), T = \\M\\ 2 = X + Y, U M = coth^(M) = ^|=. 

Since sin$M = ^= = and cos 9m = Jj^ X 4 > * ne inequalities in (7.6) can be described as 
1 |sin(g M -2g)| < r 2 < 2 

£ C/ M + cos(0 M - 20) " Q 2 " (Z7 M + cos(0 M - 20)) v / T 2_r 4' 1 ' J 

Denoting 5m = ^ — 0, from the first and last fraction in (7.7) we infer | sin 25m I *C i- Therefore 
5m is close to 0, or to ±|. When 8m is close to we have | tan 5m I *C |5m| <C | sin 25m I <C <f • 
When 5m is close to ±| we similarly have \5m T f I ^ t> which we claim is impossible. From 

I tan 5m I I sin 25m I < £ 



1 + if M ^l U M + co s 25 M ^ v 3 

l+cos 26m 

it suffices to bound from above 1H ^-os 2<5 M ' wmcri would imply | tan 5m I *C £, contradicting ^MTf I *S 
y . Since Z is a positive integer we have sin 0m 3> y- Since cos 9, sin0 > and 9m G (0, 7r), at least 
one of cos 25m +cos 9m or cos 25m — cos 9m must ne nonnegative, showing that cos 25m ^ — I cos 9m\, 
and so we have 

1 + cos 25m ^ 1 - I cos 9 m \ = 1 - yl - sin 2 9 M > ^2 • 
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As Um — 1 *C Yi i it follows that 1 _ ) y c * J s2 ^ <C 1, proving the claim. 

We have thus shown that \6m\ ^ I tan 5m I *C jr. 
Case I. Y > X. Then < e -f < f and Z = V^F - 1 < Y. Since |5 M | < y < <2 _2c °> ^ follows 
that < < | for large Q. Employing the formula tan 



— we now have 

Y —ex 



AC + BD 



C 2 +D 2 -e T q 
z 



1 + tan 9 tan 



6m 



1 



| tan 8m I 

Combining (|7.8|) with < Y -e T ~ Y ^ T ana - w ^h | y — | ^ c 2 +D^ ^ T ' we arr i ve at 

A + B g' 



C + D 9 



<^<Q- 2c °. 



(7.8) 



(7.9) 



If nonzero, the left-hand side in must be > g(c + g) ■ But g(C + J3) < q^/2(C 2 + D 2 ) < Q>/2, 
and so Cr 2c ° <C Q, contradiction. It remains that q = C + D and q' = A + B, which again is not 
possible because Q 2c ° (C + D) 2 = g(C + D) < Q\/2. 

Case II. X > Y. Then f < < f and Y < - 1 = Z. As |5 r | < Q~ 2c °, we must have 

< 5 — < ? for large values of Q. This time we have 



Y-e T 



Z 



tan 



7T _8m\ 

2 2 J 



— tan 



7T 



tan^Ml 



1 + tan 



2 

2 2 



tanf|-0^ ^ (1 + V3)|tan<5M| « ^, 



which leads (use D ^ C B ^ A) to 



C + D _ q 
A + B q' 



1 er 



Y _ C + D 
Z ~ A + B 



1 \D~C\ 
< f + (yl + B)(AC + BD) 



1 

< ~ + 



1 



(7.10) 



T (A + B) 2 T JC T 



As in Case I this leads to a contradiction because q'(A + B) ^ q'y/2X ^ Q\/2 and (A + B) 2 ^ 
Q 2c °. □ 



In summary, the nonzero contributions to (7.1) can only come from the matrices M E (3 with 
max{ Xm,Yjm} ^ Q 2c °. Next we apply Lemma [7] with r = 1 to the set 

QV 



w 2 + q 2 



O = ^M, g ',« = |(«, w) e [0, Q] x [0, g'] : \E M (q', u)| < ~, g /2 X M + g 2 >M + 2gg'Z M 

assuming max{X Af , Y M } < Q 2c °- We have A(O) < Q 1 ' ^', £(dQ) < Q^+g 7 < Or 1_c °, L = g' 5 / 6 , 
getting £n,L,q' "C Q 1- ^'' -1 / 6 and therefore 



E 



^1 

„/2 



'A(Om,^) + £ ( £ g^-Og'-Ve+A 



(7.11) 



A choice of max{A, B} and max{C, D} determines M uniquely, the total contribution of error 
terms from (7.11) to Bg(£) is 

E E E Q 1_C0 g /_1/6+e < Q ^ ff) 5/6+£ «Q (11+Co)/6+£ . (7.12) 



C^Q c o A^Q c o q'^Q/A 



A^Q c o 
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To estimate the contribution of the main terms we apply Mobius summation to the function hi (q') = 
-fA{p.M,q',i) with ||/ii||oo ^ Q 1_c ° and the change of variables (q',u,v) = (Qx,Qy,Qxz), (x,y,z) G 
[0, l] 3 , {7T§ and to find 

Vol(5^) + O c (Q 1+Co In 2 Q) + O e (Q 1+Co+£ + Q( n + c o)/6+ £ ) f 



C(2) 



(7.13) 



A A f,yAfsSQ 2c ° 

where S'jv/,^ is the three dimensional body 



1 



Sm,z ■= Ux,y,z) G [0,1] 3 : |H M (x,y)| < £, x 2 X M + y 2 Y M + 2xyZ M ^ ' (7 - U) 

With 77 = (jj) notice the equality E V Mr/(y, x) = —EM{x,y). Therefore the reflection (x,y,z) 1— > 
(y,x,z) maps Sm,£ bijectively onto S^Mr^- 

Next we give upper bounds for Vo1(»Sm,§)- Take (x,y,z) = (r cost, r sin t, z) G Sm£- The proof 
of ( 7.9 ) and ( 7.10 ) does not use the integrality of q' and q, so denoting com = xjr§ ^ < ^ M an< ^ 
W M = c+D ^ ^ M < we ^ nc ^ 1 1 ~~ T * n ^ e f° rmer case an d respectively 1 1 — <C ^ 

in the latter case. Writing the area in polar coordinates we find r 2 ^ ^ and 

Vol(5 M , 5 ) < G [0,1] 2 : 3z G [0,1], (x,y,z) G 5 M ,J) 

ar^M-^ - 2e (7 - 15) 



The bound fl7.15| ) yields 

2 Vol(5 M ,e) < 00 and £ Vol(S M , f ) < S Q~ 2c °. (7.16) 

A/e6 Me6 

A 2 +S 2 ^Q 2c o 

Prom gl3| ), ( [7161 ) and c G (§, l) we infer 

^q(0 = ^) E Vo1 (5m, 5 ) + £ (Q( 11+c °)/ 6+£ ). (7.17) 



Me6 

The volume of £m,£ can be evaluated in closed form using the substitution z = tan t: 

Vol(5 M ,f)=/ B M (Z,t)-^-, (7.18) 

where .E?m(£>£) is the area of the region 

. . .. o 1 I sin(20 - 6m)\ 9 1 2cos 2 t ] , 
rcos0,rsin# G 0, 1 2 : - • — - ^ — '} . < r 2 < - • — — L 7.19 

with 9m G (O, f ) having sin 0m = ^ 2 ^ 2 M | and C/y = ^§ 4 (for brevity we write T = Tm)- 

The following elementary fact will be useful to prove the differentiability of the volumes as 
functions of £. 
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Lemma 10. Assume G,H : K 
x+ = max{x, 0}. The formula 



v(0 



defines a C 1 map on M and 



are continuous functions on a compact set Kcl' and denote 

[ (Z-G(v)) + H(v)dv, £eR, 
Jk 

= ! II(r)rlr. 



G<5 



Using equation (7.18) we find 

2 



dt 



o 



t/2 (2/^t 2 ~^1 
de— — 



sin(2# -fl M )|/(£cos 2 i)) 



o U T + cos(2# - 6» M ) 

and applying Lemma 10 we obtain: 

Corollary 11. The function £ i— )■ Vol(jSijy ^) is C 1 . 

For a smaller range for £ we have the following explicit formula. 

Lemma 12. Suppose that £ ^ Zjy/. XTie volume of Sm,£ only depends on £ and T = ||M|| 2 : 



(7.20) 



Vol(5 M ,f) 



where A 



tt/4 



tan 



! /v 7 ^- VA - 4£ 2 COS 4 t 
2a£ cos 2 1 



4 cmc? a 



r+v / T 2 -4 
2 



2£cos 2 i V 2a 1 



Proof. The two polar curves in ( |7.19 ) intersect for | sin(2# — 6m)\ = ~7pf 4 cos 2 i that is for 0-t = 

^ ± a, with a = a(£,i) G (0, f ) such that sin 2a = -^f 4 cos 2 1. Since sin#M = 4 , the 

assumption £ ^ Z ensures that a < 6m- Thus #± 6 [0, § ), and a change of variables 9 
yields 

1 f a ( 2cos 2 t 1 | sin(2u)| 



7 -f + u 



VT 2 -4 tfr + cos(2u) £(f7 T + cos(2u)) 



du. 



The integrand is even and both integrals can be computed exactly, yielding the formula above. □ 



In particular Lemma 12 yields Yo\{Sm,{) <C providing an alternative proof for (7.15). 



7.2. Exterior arcs. Referring to the notation of Sec. 5.2, we first replace the inequalities p 2 + 
p' 2 + q 2 + q' 2 < Q 2 and p\ + gf + (i^p^ - Pf-i) 2 + (-^ft - qi-i) 2 < Q 2 in ( |5.6[ ) by simpler ones. 
Using p'q — pq' = 1 we can replace p by in the former, while pi-\ can be replaced by 



the latter. As a result these two inequalities can be substituted in (5.6) by 

„2> 



111 



1 + 



J2\ 



s$ Q 2 [l + O 



1 



1 + ^ ) («f + 



^ Q 2 1 + O 



(7.21) 



Since | 
written as 



f 7 + Oil)) anc ^ ^ (-^^ ~~ <Z£-i) 2 ^ 2Q 2 , the second inequality in (7.21) can be also 



i + ^)fe 2 + (^-^i) 2 )^Q 2 (i + o@), 
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leading to 

R Q n (£) = E E • A/ Q+0(Q 1 /2), g /,i^(£)> 
te[0,C) <z'«2 

where A/q ^ denotes the number of integer lattice points (p 1 , g) such that 

O^p'^g', O^g^Q, p'g = 1 (mod <?'), 0<K^-%_i^Q 

^Wq)^^ P q "maxigH^^ + ^-ft-i) 2 }' 
Applying Lemma [7] to the set f2 = ^^xg °^ elements (u, v) for which 



u€ [0,Q], v€ [0,<f], U 



> 0, i = 0,1, 



0<kl '[q-q 

2 i /2 ^ 



Q Q 



,2„/2 



Q 2 9 



max |«a + g* Q 2 L 2 (^, g) + Q 2 (iO^, g) - L^^, g 
with A(J2) < Qg', £(5£2) < Q, L = g' 5 / 6 , we find 



Nq2',1,k(£) 



+ O e (Qq'- 1 ^). 



This leads in turn to 
where 



^ n ^)=- M Q n (^) + ^(Q 11/16+£ )' 



^3 n (o = E E ^ 



£e[o,0 g'^Q 



For fixed integers G [1,£)j ^ e IP>£) consider the subset Tk,1£ of [0, l] 3 defined as 
< L m (x,y) = KL e (x,y) - L t _i{x,y) < l,T^(x,y) ^ O 



(s,j/,z) G [0, l] 3 : 



max 



{x 2 + y 2 ,Lj(x,y) + Lj +1 (x,y)} ^ - 



(7.22) 



(7.23) 



with Li and as in (5.4) and (5.5). 

Mobius summation applied to h^{q') = ^Aip^,^ K ^) with 1 1 1 1 00 ^ Q, and the change of variable 
(q',u,v) = (Qx,Qy,Qxz), (x,y,z) G [0, l] 3 , yields 

M ™® = 7?2) £ ( [ 9 %MMVZk*)+0(Q))=^ E Vol(T^,,), 



and so 



^(0 = ^ E VoKT^^ + O^^ 



C(2) 



(7.24) 



fc[0,£) 
*e[U) 
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To show that £ i— > Vo\(Tk^) is C 1 on [l,oo), we change variables (x,y,z) = (cos 9, sin 9, tan t) 
to obtain 

r?r/4 . , dt 



Vol(T KM )= A K 4Z,t)^-, (7.25) 
Jo cos 

where A/^(£, i) is the area of the region defined by, setting e %e = (cos 9, sin 9): 

(r^[0,lf:%(e tf )>0, %^ ^ ^ ^ 



6 max{l,Lf( e *)+Lf +1 ( e »)} 



with the notations of (5.4) and with 

L i+1 (e w ) 



Li_i(e«)Li(e«) L^){h\{^) + Lf +1 (e*)) ' 



Now notice that Ki(x,y) ^ £ when 1 ^ i ^ £, as a result of (omitting the arguments of the 
functions) 

Similarly, K\ = L ~j^ Ll ^ + L ^ Ll < ^£,K ^ £■ Thus the projection of on the first two 

coordinates is included into the union of disjoint cylinders 7k := Tk x C\T~ l Tk 2 H . . . DT - X Tk t with 
Tk = {{x,y) : K\(x, y) = k} and k = (k±, . . . , kg) G [1, £) . On each set 7k all maps Li, . . . , Li, L? + i 
are linear, say Li(x, y) = A{X + Biy where the integers Aj, -Bj depend only on . . . , ki for i ^ ^ 
and Ai + i,Bi + \ depend only on k and K. Therefore the function Fki{9) is continuous on each 
region 7k, and applying Lemma 10 we conclude that the function £ \-t VoI(Tk,i,£) is C 1 on [l,oo], 
being a sum of volumes, each of which C 1 as functions of £. 

Remark 13. The region Tk,1£ can be simplified further. For each integer J £ [1,£), the map 

: ' ^ (JL e (u,v) - Lt-i(u,v),Le(u,v)) 



is an area preserving injection on 7", since it is the composition of T in (5.4) followed by the linear 
transformation (u,v) h-> (Ju — it, t>). Note that under this map (omitting the arguments (u,v) of 
the functions below): 



Li 



1 + JLi — L^_! 



(JLe - Li-x) = 



(using Li-i + Li > 1), and by induction it follows similarly that Li — > Li-i for ^ i ^ I. Also we 
have that \Pj(it, u) = (x, y) € [0, l] 2 if and only if x = JL e - L^_ x G [0, 1] and J = . 

Let us decompose the region Tjc^ into a disjoint union of regions T^j-e.^, 1 ^ J < £, obtained 
by adding the condition [^jp] = 7- By the discussion of the previous paragraph, the map (^j, ld z ) 
is a volume preserving bijection taking Ujc j-e,(, onto Txj-£^, where 

J- s + y>l, JL e -L e _x>0, KLo-LtX), T tfK>J ^S,} 

UK ^ ■= £ [ °' : L 2 + (i^Lo - Ll f < ^ L 2 + (JL, - L,_ l} 2 < ■ 

Here L< = Lfay) and T w (x,y) = ^^Sl^zL) 2 ) + lOT + Lo^+K-L^) ' 

For a ^ 1, the transformation (\t Q , Id z ) maps bijectively the part of Uk,J;££ for which = 
a onto the part of for which \}^\ = or. Therefore Vo1(J7r- : j ; ^) = Vo^fJj^^^) and the 
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sum of volumes appearing in ( |7,26[ ) can be written more symmetrically: 

J2 Vol(T^)= Yl Vol(U KiJM ). 
Ke[i,0 K,Je[W 

As an example of using this formula, if 1 < £ ^ 2 and £ = 1, we can only have K = J = 1 and the 
inequalities JL\ — Lq > 0, KLq — L\ > cannot be both satisfied, so E/i l-i^ is empty. Therefore 
the only contribution from the T bodies in (7.26) comes from 2i o ^ if £ € (1, 2]. 

We can now prove the main theorem regarding the pair correlation of the quantities {tan(# 7 /2)}. 

Theorem 2. The pair correlation measure i?f exists on [0, oo). It is given by the C l function 



Rn 



3C(2) 



\Afes te[o,£) Ke[i,() / 



(7.26) 



where the three-dimensional bodies Sm.£ °> t £ defined in (7.14) and the bodies Tk/£ are defined in 
(733}. 

Proof. If co 6 (5, l) and G(£) denoting the sum of all volumes in ( |7.26[ ), by (7.17), (7.24), we infer 

Q 2 



C(2) 



G(£) + Ote{Q 



(ll+co)/6+e^ 



(7.27) 



The function G is C 1 on [0,oo) as a result of £ 1— > VoI(Sm,() being C 1 on [0,oo), and of £ 1— > 
Vol(T KM ) being C 1 on [l,oo). Corollary [Z] and §L27\ ) now yield, for /3 e (§, l), 



C(2) 



G(C + 0(Q 2 " 3/3 )) +G(0(Q 2 " 3/3 )) + C\ £ (Q 1+/3 lnQ + Q 



)(ll+co)/6+e 



Employing again the differentiability of G, as well as G(0) = 0, and taking /3 



leads to 



co 



C(2) 



G(0 + o?, £ (Q 23/ 



Equality (7.26) now follows from ( |7.28 ) and Corollary [8j 

8. Pair correlation of {# 7 } 
In this section we pass to the pair correlation of the angles {9-y}, estimating 

:= #{(7,7) £^:0<e y -(9 7 <^J, 
thus proving Theorem [J] In §8.2| we also deduce Corollary [TJ 
8.1. Proof of Theorem [TJ Define the pair correlation kernel F(£,t) as follows 



£, this 

(7.28) 
□ 



A/e6 



*"6[U) 



where I?Af(£,i), are the areas from (7.18), (7.25), so that by (7.28) we have 



C(2) 



cos 2 (i) 



(ll+c )/6+e 
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Proposition 14. We have the estimate 

Q 2 



(8.2) 



Before giving the proof, note that Theorem [T] follows from the proposition as Q — > oo, taking 
into account the different normalization in the definition of TZq(^), 72q(£), and defining, in view of 
dO and dO): 



Bm(0 



Bm 



2 cos 2 t ' / cos 2 i 



.4 



Kj 



2 cos 2 t ' / cos 2 t 



•3) 



From the definitions of 7?&f(£,7), Ax/(C,t) in the equations following (7.18), (7.25), it is clear that 
B M{ 2c l s i v t) = 5 M (§,0) cos 2 t, A Kfi ( 2 J sit ,t) = A Kj0 (^,0) cos 2 t, hence one has 



7T 



7T 



(8.4) 



which together with ( 7.20[ ) yields the formula for Bm{C) given in Theorem [TJ Note that the range 
of summation in Theorem [T] restricts to K < £/2, I < £/2, compared to the range in (8.1). Indeed, 
from the description of Ag^ ( 2 a . , t) following (7.25) we see that I < Ti x *s £/2, while for K we 



have K < 



+ 



KL t -K t _ x 



< ^e,K § > an< ^ similarly for t = 0. 



Proof. Consider 7 = [a,/3) with N = [Q d ], \I\ = N' 1 ~ 7+ = [a - Q~ d ',/3 + Q~ d '], 7~ = 

[a + Q~ d , /3 — Q~ d ] where 0<d=^<d' = j^<l. Partition the interval [0, 1) into the union of 
N intervals Ij = [ay, ct/+i) with \Ij\ = N as above. Associate the intervals 7- to 7j as described 
above. Denote 



:={(7,7')eiK 2 Q:777}> 



7 ^ Q 2 ' 



*( 7 ),*( 7 ')e7 



< ^Q(e) := #{(7,7) € ^ : < fly - 7 < J,, *( 7 ) € /} 
tt*Q(0 := #{(7,7) G ^ : CK - tf( 7 ) < ^, ^(7)^(7) G / 

: = #{(7, 7') G 9^ : < *(-/) - *( 7 ) ^ ^, 7-, 7+ G 7 
== #{(7, 7) 6^:0^ *(-/) - *(7) < 



i 

Q 



2 , 7- , 7+ G 7 



Expressing fly — fl 7 and ^(7') — ^(7) by the Mean Value Theorem we find 



(l + /? 2 )g 
2 



• 5) 



Lemma 15. The following estimates hold: 

N N N 

(») ^ Q (0=<i(0 + O(Q 1+d 'ln 2 Q). 
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IS 



Proof. The first inequality in (i) is trivial. For the second one note first that the total num- 
ber of pairs (7,7') with < 9y — # 7 ^ £Q~ 2 and qq' ^ Q d ' , with 7_ = |, 7+ = 
< f Q d (Q d ' lnQ)(Qln<2). For 7 with gg 7 > Q~ d ' use *( 7 ') - < tf(V) - *( 7 ) s: +7 ^ 
so ^(7') G /j". The proof of (ii) is analogous. 



Returning to the proof of Proposition 14 Lemma |15| and (8.5) yield 

N 



1 + W +Ce(QU / W 



□ 



(8.6) 



To estimate 7£jq(£) we repeat the previous arguments for a short interval I as above. Adding 
everywhere the condition 7 _, 7+ G / we modify TZq by ^7 q and Rq by i?j q in Lemma [5j T^-q by 
7^.jq and -Rg n by -Rpg in Lemma |6j The additional condition ^ G / is inserted in (|7.2|). The 



condition ^ p' ^ (/ is replaced by q'a ^ p' < q' f3 in (7.4), and ( 7.22[ ), and ^ p ^ g is replaced 
by c/a ^ p < qj3 in (7.4). The condition v G [0,q'] is replaced by v G [q'a,q'(3) in the definition of 
^M,<?',£) and The bodies Sm,£ and Tk/£ are substituted respectively by Si t M£ and T^k,1£ 

after replacing the condition 2 G [0, 1] in their definition by z G [a,/?). The analogs of (7.18), and 
( 735) hold: 

Vol(^, M , f ) = / 2-, Vol(T I>KM ) = / Aj^(£,t)— ^ 

7 1 cos z t J I cos z t 



(8.7) 



The approach from Section [7] under the changes specified in the previous paragraph leads to 



^(0 = ^(0 + ^(0 



C(2) 



^ , j. s dt 

-1 r COS^ i 



+ ? ,e(Q 23/12+£ ), 



with the pair correlation kernel F(£, t) defined by ( |8.1[ ). We also have 

r% q (o = (o + o^(Q 23/12+£ + Q 2 - d '; 



The analogs of Lemmas |5j |6| yield upon (8.8) and ( |8.9 
Q 2 



<'q(0 



C(2) 



/ ^ + o(Q- 1/3 M)^ + o c , £ (Q 23 / 12+£ ) =K*f Q (0. 

Aan- 1 / 'cOS^t 

The analog of Corollary [4] and (8.10) yield 

= n%{t + o{Q-y±)) + k^(o(q-^)) + o(q^) 



(8.9) 
(8.10) 



C(2) 



/ (F^ + OiQ-^^+FiQ-^^))) 

J tan- 1 1 V 7 



(ft 

cos 2 i 



.11) 



K 



/+,< 



>(0- 



As shown in Section [7] the function F is C 1 in £, thus ( 8.11[ ) gives actual!^] 



<£(0 ^ 



C(2) 



/ x mt)^ t+ O,AQ 23/12+£ )=^\ Q (0. (8-12) 

•/tan" 1 1 Lob 1 



Lemma 15 (i), dglgh, and F G C^O.oo) yield 



>(0- 



.13) 



^The argument from Section [7] applies before integrating with respect to t on [0, j] , showing that F is C 1 . 
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Let also ujj = tan -1 ay. From (8.13) and (8.5) we further infer 

Q 2 m+i /(l + a 2)^ \ dt 
pi 



23 



C(2) 



Q 2 P+ 1 ^/(l + a 2 +i)^ A 



C(2) 



(8.14) 



* ^7 + ^ £ (Q 23 / 12+£ + Q 2 - rf '). 
' cos z t 



Employing also 



F 



cos 2 i 



Wj+l 



F 



(1 + tan 2 f)f 



,t 1 + 0(wj+i - Wj) 



cos 2 t 



and (wj+i — Wj) < Q we find 

q2 r u j+1 



C(2) 



1 / (1 + tan 2 i)£ \ 



' t )c^ + ^^ 23/12+£ ) = ^;- Q ^ ) - (8 - 15) 



Finally Lemma uM{i) and (8.15) yield (8.2). 



□ 



8.2. Explicit formula for gf. Next we compute the derivatives B' M (£), thus proving Corollary 
[l] We also obtain the explicit formula (8.16) for on a larger range than in Corollary [TJ after 
computing the derivative A' KQ (^). 



Lemma 16. For M G (3, let T = Tm, Z = Zm as in (3.1). The derivative B' M (£) is given by: 



71 ]n( T+ ^ TT ~- 1 
4£ 2 n VT + \JT 2 - 4 - £ 2 



In 



(T + a/T 2 ~^4) 2 (T - ^T 2 - 4 - £ 2 
(4 + 4Z 2 )(T + \JT 2 - 4 - 



8£ 2 

vr , /(r+v 7 ^ 2 ^) 2 

,,,( 4 + 4Z 2 



ift^2Z 

if2Z ^ £ ^ VT 2 -4 
> VT 2 - 4. 



Proof. Using (8.4), we proceed as in the proof of Lemma 12 

rir/2 



TT 



l 



sin(20-0 M )| 



de, 



^W^l U T + cos(2# - M ) U T + cos(26 - 9 M ) , 

we obtain: 



where Ut = 4 and 6* m £ (0, f ) has sin#jv/ = /^f ^ ■ Applying Lemma 



k_ f |sin(2fl-fl M )| 
4£ 2 J i U T + cos(20 - 9 M ) 



10 



de, 



with J = he (0, f) : |sin(20 



Ml < Tffcl}- Clearly / = (0, f ) when £ > VT 2 - 4, and if 
f < VT 2 - 4, let a = a(0 G (0, |) such that sin 2a = ^^j- Then 



^ < 2Z a < M /2 = [0 M /2 - a, M /2 + a], 

2Z < £ ^ V^ 2 -4 a G [6» M /2, tt/4] I = [0, A/ /2 + a] U [vr/2 + M /2 - a, tt/2], 

■11 

.01 



and the integral is easy to compute. The last case appears in the figure below for M = (}>]) and 



£ = 2.lG (2,^5). 



□ 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 6. The area B M (£/2, 0) when 2F ^ f ^ yjT 2 - 4. 



A similar computation using (8.4) gives the formula: 




TT 
4^2 



ln(l + K 2 ) + In 
ln(l + i^ 2 ) 



(l + x 2 )(l + (x 2 -iQ 2 ) 
(1 + x 2 )(l + (xi - K) 2 ) 



if £ ^ 2K 

if £ £ [2K, KVK 2 + 4] 
if £ ^ £Vif 2 + 4, 



where x 2 > x\ are the roots of x 2 {^ + 2K) - 2xK(£ + K) + £{K 2 + 1) - 2K = 0. By the last 
paragraph in Remark 13 the body ?i i £ is empty, so Ai i(£) = 0, and we have an explicit formula 
on a larger range than in the introduction: 



9l 



3 

4tt 



< i < 4. 



(8.16) 



We can now explain the presence of the spikes in the graph of in Figure [TJ The function 
B' M (£) is not differentiable at £ = 2.F and y/T 2 — 4, while the function A^o(£) ^ s n °t differentiable 
at £ = 2ivT and (X 2 + 2) 2 — 4. At the point £ = \/5, two of the functions B' M (£), as well as 
A'i (£), have infinite slopes on the left, which gives the spike on the graph of g^{x) at x = /-\/5. 
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